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Abstract 

We study the pair correlation function for a variety of completely integrable quantum maps in one 
degree of freedom. For simplicity we assume that the classical phase space M is the Riemann sphere 
(DP^ and that the classical map is a fixed-time map exptEn of a Hamilton flow. The quantization 
is then a unitary N x N matrix Ut,N and its pair correlation measure Pjl' gives the distribution of 
spacings between eigenvalues in an interval of length comparable to the mean level spacing (~ ^/N). 
The physicists' conjecture (Berry- Tabor conjecture) is that as ~+ oo, pj^^ should converge to the 
pair correlation function pPOissoN = ^ ^ Qf ^ Poisson process. For any 2-parameter family of 
Hamiltonians of the form Ha^/s = ot(j){I) -I- 131 with 0" 7^ we prove that this conjecture is correct for 
almost all (a, 13) along the subsequence of Planck constants Nm. ~ [m(log m)*]. In the addendum to this 
paper [Z. Addendum], we further show that for polynomial phases cj> the a.e. convergence to Poisson 
holds along the full sequence of Planck constants for the Cesaro means of p2,t,a,p- 
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Introduction 



In this paper we shall be concerned with the fine structure of the spectra of some completely integrable 
quantum maps in genus zero, that is, with quantizations of integrable symplcctic maps x on the Riemann 
sphere M = CP^, equipped its standard (Fubini-Study) form lo of integral area. For any positive integer 
N, {<CP^, Nlu) is quantized by the Hilbert space Hn = T{L"^) of holomorphic sections of the Nth power of 
the hyperplane section line bundle. The quantum system then consists of a sequence of unitary operators 
{U-^^n} on Hn, a Hilbert space of dimension N. For simplicity, wc restrict attention to quantizations Ut^N 
of Hamilton flows xt = expfSjj where the Hamiltonian H has no separatrix levels. Our interest is in the 
semiclassical asymptotics {N oo) of the pair correlation function p^,^-* and number variance of 
the quantum systems. We first show that the time-averages of these objects tend to the Poisson limits 

as N 00. This is consistent with the Berry- Tabor conjecture [B.T] that eigenvalues of completely integrable 
quantum systems behave like random numbers (waiting times of a Poisson process). However, it is only a 
weak test of the conjecture since the averaging process itself induces a good deal of the randomness. A 
much stronger test is whether the variance tends to zero. For special 2-parameter families of Hamiltonians 
Ha,p = a(f)(I) + /?/ (sec §2 for the definition) wc show that the variance tends to zero at a power law rate. 
This implies that the individual systems are almost always Poisson along a slightly sparse subsequence of 
Planck constants. In an addendum [Z. Addendum] we will further show that when ^ is a polynomial, then 
the Cesaro means in A'' of the pair correlation function p2^t\,0 ^^^d almost always to p2^^^^^^ . 

Before describing the models and results more precisely, let us recall what the level spacings problems 
are about. In the quantization H H^-^^ of Hamiltonians on compact phase spaces M of dimension 2/, 
the 'Planck constant' is constrained to the values h = 1/N and the spectrum of H^^^ consists of djv ^ 
eigenvalues {Xnj} in a bounded interval [minH,maxH]. Similarly, the spectrum of a quantum map U-^,n 
consists of djv eigenvalues {e'^"'^} on the unit circle S^. The density of states in degree N 

has a well defined weak limit as A^ oo which may be calculated by standard methods of microlocal analysis 
(§2). According to the physicists, there also exist asymptotic patterns in the spectra on the much smaller 
length scale of the mean level spacing -j^ between consecutive eigenvalues. The pair correlation function p2, 
for instance, is the limit distribution of spacings between all pairs of normalized eigenvalues cInXnj- The 
length scale ^ is usually below the resolving power of micrlocal methods. Hence; the problem of rigorously 
determining the limit , or even of determining whether it exists, has remained open for almost all quantum 
systems. The sole exceptions are the cases of almost all flat 2-tori [Sa.2] (see also [Bl.L]) and ZoU surfaces 
[U.Z]. For other rigorous results on level spacings for Laplacians on surfaces with completely integrable 
geodesic flow, see [S][K.M.S][B1.K.S]. 

On the other hand, there exist numerous computer studies of eigenvalue spacings in the physics literature 
which indicate that limit PCFs often exist. The following conjectures give a rough guideline towards the 
expected shape of the level spacings statistics: 

• When the classical system is generic chaotic, p2 = P2^*^^ where p2^^ is the limit expected PCF for 
NxN random matrices in the Gaussian orthogonal ensemble; 

• When the classical system is generic completely integrable, p2 = p^O/SSOJV ;= i _|_ That is, at least 
on the level of the PCF, the normalized spacings between eigenvalues behave like waiting times of a Poisson 
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process. The term Sq comes from the diagonal, while the term 1 reflects that any spacing between distinct 
pairs is as likely as any other. 

These conjectures should not be taken too literally, and indeed cannot be since the term 'generic' is 
not precisely defined. Our main purpose in this article is to test the Poisson conjecture against quantized 
Hamilton flows in one degree of freedom on the compact Kahler phase space CP^. Of course, they are 
necessarily completely intcgrable. It might also be suspected that quantized Hamilton flows in one degree of 
freedom are necessarily trivial, but this is not the case: as will be seen, PCF's of toral completely integrable 
systems on CP^ are almost always Poisson along a slighty sparse subsequence of Planck constants. It should 
also be recalled that many of the model quantum chaotic systems, such as kicked tops and rotors and cat 
maps, take place in one degree of freedom and still defy rigorous analysis [Iz][Kea]. 

The quantum maps studied in this paper thus join a growing list of integrable quantum systems whose 
level spacings have been shown rigorously to exhibit some degree of Poisson statistics. On the other hand, 
it is clear that not all of the quantum maps in our 2-parameter families exhibit Poisson behaviour (see the 
Appendix for a counterexample). Rather, our results tend to corroborate the (still rather vague) picture 
that (some level of) Poisson statistics occurs almost everywhere in an n-parameter family of non-degenerate 
integrable systems, but that a dense exceptional set of non- Poisson systems occurs as well. This probabilistic 
revision of the Berry- Tabor conjecture seems to have been first proposed by Sinai in his study of a closely 
related lattice point problem [S]. It was also stated clearly by Sarnak [Sa.2] in his proof that almost every 
flat 2-torus has a Poisson PCF but that a dense residual set of flat tori had no well-defined PCF. 

We should also emphasize that all of the rigorous results at the present time on eigenvalue spacings of 
completely integrable systems only pertain to the 2-level correlation function; neither the k- level correlation 
functions for fc > 3 nor the (nearest-neighbor) level spacings distribution have been proved to be Poisson. In 
this connection it is interesting to recall a suggestion of Sinai concerning the 'degree' of Poisson behaviour 
of typical members in an n-paramctcr family of quantum systems (e.g. as measured by the largest k so that 
the A:-level correlation function is Poisson) : namely, that the typical degree could depend on the number 
n of parameters in the system. Thus our 2-parameter families (and the 2-parameter family of flat 2-tori) 
exhibit Poisson PCF's a.e., but it is unknown whether their higher-level correlation functions or level spacings 
distribution are also Poisson. 

Now let us be more precise about the models we will study. In the usual Kahler quantization of (CP^, iVw), 
Hn may be identified with the space Vn of homogeneous holomorphic polynomials f{zi, Z2) of degree N on 
C^. A classical Hamiltonian H e C°°(CP^) is then quantized as a self-adjoint Toeplitz operator 

:= UnHUn -.Hn^Hn, Vjv,i ^ TInH^j^j 

where Hjv is the (Cauchy-Szego) orthogonal projection on Hn (§!)• Hence the quantum Hamiltonian system 
amounts to the eigenvalue problem: 

H''^^(j)pfj = Xnj'Pnj, — ||-ff||oo < ^N,l < Xn,2 < • • • Xn,N < ll^^lloo- 

For a fixed value h = of the Planck constant, the distribution of normalized spacings between all 
possible pairs of eigenvalues of H^^^ is given by the Nth pair correlation 'function' (measure) 

1 ^ 

Here, the eigenvalues are rescaled, Ajv,j ^^N,j to have unit mean level spacing, i.e. so that N{XN^i+i — 
Ajv,i) ~ 1 on average 
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Our first result gives an explicit formula for the limit pair correlation function 

of a quantized Hamiltonian H^^^ . It is of a similar nature to the pair correlation function for a ZoU Laplacian 
([U.Z]) and involves dynamical invariants of the classical Hamiltonian flow expiSij generated by H on the 
classical phase space. Under some generic hypotheses (which will be stated precisely in §2), the formula is 
given by: 

Theorem A For the generic H e C°°{M), the limit pair correlation function for the system H^^^ is given 
by: 

P2{f) = Vm + E E E / f{kT^{E))T^{EfdE 

where: 

(i) V = vol{{zi, z2)eMxM: H{zi) = H{z2)}; 

(ii) {c^} is the set of critical values of H ; 

(Hi) For a regular value E e (ci/, c^+i), H^^{E) is a union of periodic orbits {'yi}JJ^i'' ofex.pt'E.H and T^{E) 
is the minimal positive period of the jth component. 

It follows that the pair correlation function of quantized Hamiltonians in one degree of freedom is quite 
deterministic. On the other hand, the eigenvalues of the associated quantized Hamiltonian flow are much 
more random. 

Before describing the results, let us recall the definition of a quantum map and of its pair correlation 
function. Suppose that Xo is a symplectic map of a compact symplectic manifold {M,u)). It is called quan- 
tizahle if it can be lifted to a contact transformation x of the prcquantum bundle tt : {X,a) {M.uj), 
where da = 7r*a;. We are mainly interested here in Hamiltonian flows xt and these are always quantizable 
(§1). We then deflne the quantization of the map Xo to be 

Ux,N •= n^vCT-^TxHiv : Ti-N — * Wat 

where is the translation operator by x on T-Lm and cr^ G C°°(M) is the 'symbol', designed to make U^^n 
unitary. All of the usual quantum maps, e.g. 'cat maps' and kicker rotors can be obtained by this method 
[Z]. 

Since the eigenvalues lie on the unit circle, the rescaling to unit mean level spacing leads to the (period) 
pair correlation functions 

dpr(^) ]v E E ^(^ - ^(^^.' - ^^.^ - ^^)) = E E ^^^^'^ - ^^.^ - 2^^^ - |f)- 

The large N behaviour of dp'^^ is quite 'random' in general because the rescaling destroys the Lagrangean 
nature of U^,n- The question is whether there is some asymptotic pattern to the randomness. 

As mentioned above, we will restrict in this paper to a simple but reasonably representative ease of the 
question, namely to Hamiltonian flows generated by perfect Morse functions H on CP^. The reason for 
restricting attention to CP^ is that it is the only symplectic surface carrying a Hamiltonian action (i.e. 
it is a 'toric variety'), namely the usual rotation of the sphere about an axis. The moment map is known as 
an action variable /. Any perfect Morse function may be written as a function H ~ (t>{I) of a global action 
variable. Any toral action can be quantized and in particular / can be quantized as an operator /(^) whose 
spectrum lies on a one dimensional 'lattice' {j^ '■ j = —N . . . N}. It follows that H is quantized as an operator 
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of the form (/>(/) and its flow can be quantized as a unitary group of the form Ut^N = ^Ne^*'^^''"^IlN, where 
Af equals to TV on Tijv. Hence the cigenanglcs have the form tN(j)[jj) and the asymptotics of the PCF can 
be reduced to the study exponential sums of the form 

N 

The Poisson conjecture is essentially that these exponential sums behave like random walks. It is too difficult 
to analyse the individual exponential sums, but we can successfully analyse some typical behaviour in families 
of such systems. The first result is about the mean behaviour as the t parameter varies. 

Theorem B (a) Suppose H : M ^ IR is a perfect Morse function on CP^. Then: the limit PCF p2,f and 
number variance I^2,t{L) for Ut,N are Poisson on average in the sense: 

lim J- [''ffdt = p^O^'':=l + 6o 

and 

lim -1- ri]W(L)di = Sr^^^(L):=L 
N—Kx — a J ^ 

for any interval [a, b] of IR. 

This result applies to the case of linear Hamiltonians and their Hamilton flows, whose pair correlation 
functions are clearly not individually Poisson (cf. §6). For Poisson level spacings, we make some further 
assumptions on the Hamiltonian (or phase (j)). Our main result concerns the mean and variance of a 2- 
parameter family of Hamiltonians: 

Theorem B (b) Let I denote an action variable on (DP"'^ and let i?a,/3 = + (31 with \6"\ > 0. 

Denote by Pj^^ q, the pair correlation measure for the quantum map U(^t^a,i3),N = ^xp{itJ\fH(^a,i3;N))- Then 
for any t ^ 0, any T > and any f G S{JR.) with f € C^{JR) we have 



Thus, the mean pair correlation function in the family is Poisson and the variance tends to zero at the 
rate • Following [Sa.2], we conclude: 

Corollary Let Nm = [rnQogm)"^]. Then, for almost all {a, (3), P2;[t,a,fi) ^ 

It would be interesting to study quantizations of Hamilton flows in the case where the Hamilton had 
saddle levels, as must happen if the genus is > 0. It would also be interesting to study completely integrable 
maps which are not Hamilton flows. We hope to extend our methods and results to these cases in the future. 

Acknowledgements This article was completed during visits to the Australian National University and 
to the Newton Insitute. In particular, we thank A.Hassell and Z.Rudnick for comments on the proof of 
Theorem B and J.Marklof for discussions of incomplete theta series. We particularly thank M.Zworski for 
suggesting a simplification of the proof of Theorem 5.1.1 
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1 Toeplitz quantization 



We now review the basics of Toeplitz quantization on CP^. For futher background on Kahler quantization, 

we refer to [G.S]: for general Toeplitz quantization we refer to [B.G][Z]. 

Toeplitz quantization is a form of Kahler quantization, that is, of quantization of symplectic manifolds 
in the presence of a holomorphic structure. The basic idea is that the quantum system is the restriction of 
the classical system to holomorphic functions. 

To be more precise, let [M, uj) be a compact Kahler manifold with integral symplectic form. Then there 
is a positive hcrmitian holomorphic line bundle L M with connection 1-form a whose curvature equals 
UJ. In Kahler quantization, the phase space (Af, uj) is (quantized as the sequence of finite dimensional Hilbert 
spaces r(L®^) where T denotes the holomorphic sections. In Toeplitz quantization, these spaces are put 
together as the Hardy space H'^{X) of CR hmctions on the unit circle bmidle X in L*. 

Thus, the setting for Toeplitz quantization is a compact contact manifold (X, a) whose contact flow 

/ : X ^ X, /*a = a (1) 

defines a free 5''^-action with quotient a Kahler manifold M whose Kahler form uj pulls back to da. The Kahler 
structure on M also induces a CR structure on X. The corresponding Hardy space H'^{X) is the space of 
boundary values of holomorphic functions on the disc bundle of L* which lie in L'^[X). The orthogonal 
(Cauchy-Szego) projector H : L'^{X) — > H'^{X) defines a Toeplitz structure on X in the sense of [B.G]. 
Prom the symplectic point of view, H'^{X) is viewed as the quantization of the symplectic cone 

S = {(a;, ra^) : r e IR+} C T*X - 0. 

To be precise, H is a Hermite Fourier integral operator with wave front set on the isotropic submanifold 
S* := {{a,— a) : cr e S} C T*(X x X). The CR structure corresponds to a positive definite Lagrangean 
sub-bundle A of TS-*-, the symplectic normal bundle of S. The vector fields generating A annihilate a ground 
state eA in the quantization of the TE-"-. The symbol of H is the orthogonal projection it = e^® e.\ onto 
this ground state. For a detailed account of these objects we refer to [B.G]. 

1.1 Toeplitz quantization in genus zero 

In the case of M = CP"'^, the contact manifold X may be identified with SU{2) and the Hardy space H'^{X) 
may be identified with the space of lowest weight vectors for the right action of SU{2) on L'^{SU{2)). To 
make the Toeplitz theory more concrete, let us recall how these identifications are made. 

We first recall [G.H, §1.3] that the holomorphic line bundles over CP^ are all powers H^^ of the hy- 
perplane bimdlc H CP^, whose fiber over V € CP""^ is the space V* of linear hmctionals on the line 
thru V. The Chern class of H is the Fubini study form cups, which generates H^{<CP^). The holomorphic 
sections Hn of H are given by the linear functional L on by setting sz,(V^) = L\v- More generally, the 
holomorphic sections of correspond to homogeneous holomorphic polynomials of degree on C^. 

The associated principal to H is evid(^ntly the unit sphere S'^ C C^ which we identify with SU{2). As 
the boundary of the unit ball B C C^, it has a natural CR structure. The associated Hardy space is the usual 
space of boundary values of holomorphic functions on B. Under the action e'^^{zi, z^) = (e'^^i, of 
CP^, it is evident that the holomorphic functions transforming by e*^^ are given by homogeneous 
holomorphic polynomials of degree A^. The Cauchy -Szego kernel is given by IIn[z,w) = {z,w)^ . 

We also recall that the irreducible representations Hn of SU{2) arc given by its actions on homogeneous 
polynomials. By the Plancherel theorem, L'^{SU{2)) = e?^^i'HAr(8)7i^. The CR structure induced on SU(2) 
by the identification = dB = SU(2) is equivalent to that given by the lowering operator L_ for the right 
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action. The Szcgo projector Iliv is then the orthogonal projection onto Hn <S) i^N where Vjv is the lowest 

weight vector in 7i^. 

Below we will often refer to an action operator /'^^ on CP^. It may be identified with the Planck constant 
times any generator (e.g. L^) of a Cartan subgroup of SU{2). Thus its eigenvalues in Hn are the weights 

J_ 

1.2 Quantum maps 

Symplectic maps Xo on CP^ may be quantized by the Toeplitz method as long as x lifts to a contact 
transformation x of {X, a). The Toeplitz quantization is almost the translation operator by x compressed 
to the Hardy space H'^{X). Since does not usually usually preserve H^{X), IIT^II is not generally unitary; 
to unitarize it one must equip it with a symbol. In [Z] it is described how to construct a symbol on M 
for any quantizable symplectic map on any compact symplectic M so that 

:= Ua^T^U 

is unitary. We will describe the symbol in some detail in §2.3. It automatically commutes with the action, 
so is the direct sum of the finite unitary operators, f/x,jv on Ojv- We define C/x>JV to be the quantization of 
Xo with semiclassical parameter 1/A''. Its eigenvalues have the form 

Sp{U^,N) = {e^^''"- :j = l,...,dN} (2) 

where = dimH^{N) = N. 

Consider now the case of Hamilton fiows Xot = exptS// on a general symplectic manifold {M,co). 

Proposition 1.2.1 Hamilton flows are always quantizable. 

Proof What needs to be proved is that exptEn always lifts to a contact flow Xt '■ X ^ X. Equivalently 

that Eh lifts to a contac;t vector field, say Xh- We prove this by lifting expfSjj to a homogeneous Hamilton 
flow expiS// on the symplectic cone S. Let us define the function 

r : S — > IR"*", r{x,rax) = r. 

Thus, T, = X X JR^, X = {r = 1} and the IR'*' action is generated by the vector TZ = r£p. 

The natural symplectic structure a; on S is the restriction of the canonical symplectic structure ojt'X on 
T*X, which is homogeneous of degree 1. Denoting hy tt : X ^ M the projection, we have: 

ui = r-K*LOM + dr f\a. (3) 

The proof is simply that (jJt*x = dar'x where ar'x is the action 1-form. This equation restricts to S 
where ax'X = ra. Taking the exterior derivative gives the formula. 

Now return to if G C°°(M) and consider the Hamiltonian H{x,r) = rix*H{x) on S. It is homogeneous 
of degree 1 so its Hamilton vector field = w~^{dH) is homogeneous of degree zero and then its Hamilton 
fiow exptSjj is homogeneous of degree one. We claim that (i) the flow preserves X; and (ii) its restriction 
Xt to X is a contact flow lifting exptEn- 

Indeed, we have 

i-s^w = d{rH) = rdH + Hdr = rt^^uiM + i-Sgdr A a 
= ri-HjjWM - a{Ejj)dr + dr{Ejj)a. 
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Since all terms except dr{^fj)a are d^-independent we must have drCBg) = 0. Here d6 denotes the vertical 
one form oi X. It is then obvious that 

-a{Eg) = H, iSfi^M = dH. 

The second equation says that Eg projects to Eh, i.e 5^^ is a lift of Eh- Since 

C~ .a = c~.da + d(i~ _ a) = dH — dH 

we also see that Eh is a contact vector field (here, L is the Lie derivative). ■ 

2 Density of States, pair correlation function and number variance 

2.1 DOS 

Before considering the pair correlation function, we first describe the limit density of states (DOS) of quantum 
Hamiltonians and quantum maps in the Tocplitz setting. They can be easily determined from the trace 
formulae of [B.G] and indeed the calculation is carried out in [Z, Theorem A]. Let us recall the results. 
In the case of Hamiltonians, the DOS in degree N is defined by 

dp'^^\\):=—Y.^{\-\^^i). (4) 

By [B.G, Theorem 13.13] we have: 
Proposition 2.1.1 The limit DOS is given by 

Poif) = f f{H)Lo if G C(IR)). 
Jm 

In the case of quantum maps the DOS in degree N is defined by 

dpW(^):=-L^5(^_e2"«"..) z€S\ (5) 

i=i 

The limit DOS (3o is determined in [Z, Theorem A] and depends on whether the classical map is periodic 
or aperiodic (i.e. the set of periodic points has measure zero). 

Proposition 2.1.2 Let x be a symplectic map of {M,ui). 

(a) In the aperiodic case, (3 = CgdO where Co is the constant {Jj^ad^) with a the symbol ofUy^^. 

(b) Ifx'^ = id, then (3 is a linear combination of delta functions at the kth roots of unity. 

2.2 The pair correlation function and number variance 

We recall here the definitions of the pair correlation function and number variance for quantum maps U-^^n 
in / degrees of freedom. Then dimHN = d^ and the spectrum has the form Sp(i7^,7v) = {e*^"^ : j = 

1, . . . , div}. The spectrum may be identified with the periodic sequence {O^j + 27rn : n E Z, j = 1, . . . , rfjv} 
and then rescaled to given a periodic sequence of period N and mean level spacing one: {djv^jvj + 27rndjv : 
n e Z,j = 1, . . . ,djv}. 
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Definition 2.2.1 The pair correlation function of level N of a quantum map in f degrees of freedom is the 
measure on IR given by 

dp2^\x) ■= -1— ^2 ^2 ^i^N{ONj - Onj) + 2TmdN - x) (6) 

j,k=l n^TZ 

The limit pair correlation function is then: 

dp2^\x) = w — lim dp2^\x). 

iV— >cx) 



We often write the integral j-^fdp^^dx as 



H/) = :r E E fi'^NiONj - 0Nj) + 2nndN). 
By tlie Poisson summation formula we have: 

P^H/) = i E /(|^) E e^^(^--^-^) = iE/(^)l™x,ivr (7) 

i^2Z j,k=l '^N ^ «W 

A closely related spectral statistic is the number variance for the quantum map U-)^^n- It is defined as 
follows (cf. [Kea]): First, define the density of the scaled eigenangles by 

pW{e) = g ^ - dNONj + 2ndNn) = E [rf^ E e^^'''--e-^^'''/''-] 
= l + J-yTrUij,e-^-'''/''-. 

Then: 

Definition 2.2.2 The number variance of U-^,n is defined by: 



1 /•« rX-t-l^/^ 

^i''\L) = - 1/ ps{y)dy-L\'dx 



We observe that Sg^^ [L) is similar to p^'^ (/) for / = except that the t = Q term has been removed. 
2.3 Asymptotics of traces and exponential sums 

Before getting down to our specific models, let us make some general remarks about the exponential sums 

S{NJ.) ■■=TrUlj,. 

First, the traces Trll^ j^ have complete asymptotic expansions as ^ oo. To state the results, we need 
some notation. Recall that the 5^ action on X is denoted cj)^. For each £ put 

= {^jmod27r : Fix(/J o x) 7^ 0}. (8) 
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Assuming (as we will) that the maps have clean fixed point sets, the set &x,^ finite and Fix{(t)'^i o x) is 
a conic submanifold of S. We denote its dimension by ej and its base Fix(^^-i o x) r\^hy SFix(0^j o x)- 
The trace asymptotics then have the form: 

Proposition 2.3.1 For each £ ^ 0, 

oo 

for certain coefficients a^j^r ■ The leading coefficients are given by: 

where dfi-^^ is the canonical density on SFix{(j)^i °X^)- 

Before sketching the proof, let us describe in more detail the ingredients that go into the principal 

coefficients. Roughly speaking, the quantized map f/^.Ar involves two pieces of data in addition to the Szcgo 
projector: the map x ^iid the symbol a^. We recall from [Z] that the scalar principal symbol of is 
given by cr^ = (x*eA,eA)~^ where e^ is the section of the 'bundle of ground states' corresponding to 11. 
See [B.G, §11] for the definition and properties of the ground states (normal Gaussians). Since x is rarely 
holomorphic, it will generally not commute with 11 and will take ca to another ground state x*c\- After 
trivializing the 'bundle of ground states' the map x* '^'•^Y be described as follows: the derivative dx defines 
a linear symplectic map dxls-^ on the symplectic normal bundle TH^ of S. The quantization of the normal 
space is a space of Schwartz functions and the quantization of dxls^ is its image M.{dx\Y,-^) under the 
metaplectic representation. Then x* = M.{dx\T,^)- 

As a Fourier integral Toeplitz operator, S I°{X x X,graph{xy), where 

graphixY = {(<T,-x(a)) : a G E}. 

Hence its symbol is a symplectic spinor on the graph. Under the natural parametrization j : E graph{x), 
the symbol may be viewed as a symplectic spinor on S and as discussed in [Z], unitarity of U-^ forces it to 
have the form 

3*(^Uy, = (x*eA,eA)"^Mo-|2 O ca O (9) 
where \da\ is the symplectic volume density on S. Hence as a symplectic spinor on graph{x) it has the form 

f^(7^ = (x*eA,eA)~^(x*|'^o-|^ o |d(7|5) O x*eA O (10) 

Now let us describe the symbolic aspects of the trace, assuming that x has clean fixed point sets. Then 

each component of _Fi.x(x|s) carries a canonical density dV and by inserting the radial vector field TZ we 
get a canonical Liouville density d/x^ = indV on the base SFix{x) of the fixed point set. Moreover, in the 
normal direction we have the symplectic spinor factor x*eA<8)e^. Taking the trace, we get the matrix element 
(x*eA, ca). This cancels the scalar principal symbol, so the symbolic trace just gives the Liouville volume of 
the fixed point set, as stated above. 

Having discussed the ingredients in the above Proposition, we now sketch the proof. For further details 
we refer to [B.G, Theorem 12.9]. 

Sketch of Proof Consider the Fourier series 

oo 

T^,,(^) = E TrUl^e^""' = TrU'^e^'^U (11) 

JV=1 
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By the composition theorem of [B.G], is a Lagrangean distribution on with singularities at the values 
9j e and with singularity degrees beginning at ^. Hence: 

oo 

'^x,i- Y'^s„i,rU'_^_^{e-ej) (12) 

where Um{6) = N"^~^e'^^^. Note that is a periodic distribution with the same singularity at ^ = 

as the homogeneous distribution (6 — 6j + iO)~~'^. 

The leading coefficients are given by the principal symbols of T^ f {0) at the singularities. From the 
symbol calculus of [B.G], the coefficients are given by the symbolic traces described above. The expansion 
stated in the Proposition then follows by matching Fourier series. ■ 

Examples Let us consider the form of the trace for quantum maps in one degree of freedom: 

(a) Suppose that xt ~ exp tEn is the fixed time map of a Hamilton flow. The fixed point set then consists 
of a finite number of level sets {H = Ej{t)}, which must be periodic orbits of period t. Pick a base point 
rrij on each orbit and lift it to a point Xj lying over rrij in X. Then the lift of {H — Ej{t)} to Xj is a curve 
which begins and ends on Tr~^{mj). The difference in the initial and terminal angle is of course given by the 
holonomy with respect to a. This holonomy angle 6j is independent of the choice of mj and of xj and Q^^t 
is the set of these holonomy angles. Then SFix Xt ° (t>^^ is two dimensional and hence e = 3. 

As will be seen below, expressing TrU^^^j^ in terms of its eigenvalues produces a classical exponential sum 
in the completely integrable case. The trace expansion produces a dual exponential sum involving dynamical 
data. The principal term can be obtained by applying the van der Corput method [G.K][H.l], i.e. Poisson 
summation followed by stationary phase. However, the existence of a complete asymptotic expansion woiild 
probably not be clear without the Toeplitz machinery. The trace expansion (or at least the principal term) 
is often used in the physics literature under the name of the Gutzwiller trace formula. Some remarks on its 
limitations are included below. 

(b) Suppose next that Xo has only isolated non-degenerate fixed point sets. Then x fixes the entire fiber 
over each fixed point. The only singularity occurs at ^ = and the dimension of Fixx equals one. Hence 
e = 2. 

Special case: Quantum cat maps Ug^N These are the most familiar examples of quantum maps, so let 
us see what the above proposition says about them. In this case, the trace can be calculated exactly and 

equals the character of the finite metaplectic representations. For the exact calculation in Toeplitz setting, 
see [Z]. For the physics style calculation, sec [Kea]. Here we do the calculation asymptotically. 

First we observe that ii g — ^ ^ is hyperbolic then it has non-degenerate fixed points at (x,^) € 

/ such that g{x,^) = {x , ^)m,odZ'^ , i.e. at the points {g — I)^^ 2Z^ . The lifted map Xg actually has no 
fixed points. However, Xg°4'^ has fixed points if and only if [{g-{x, ^), e2'^*(*+^))] = [(x, ^, 9)], where the bracket 
denotes the equivalence in the quotient space. Since g-{x, ^) = (x, ^)+(m, n) for some (m, n) G we get that 
[(a;,0 + (m,n),e2'^'(*+^))] = [{x,^,e)]. But [(a;,0 + (m,n),e2"'(*+^))] = [(x,^), e2^'(*+^)e*'^'"V^«^'«)'('"'"))]. 
It follows that ^ 

Omn = --j^{mn + uj{{x,^), (m,n)). 

Hence ^ 

6g,i = {--(mn + uj{{x, C), (m, n)) : g{x, = {x, + {m, n)}. 

The fixed point set of Xg is therefore clean and has dimension one (the fibers over the fixed points of g 
on M'^/Z'^). The fixed points are non-degenerate in the directions normal to the fibers, so the canonical 
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density is given by , ^ d9 where d6 is the invariant measure on the fibers. Furthermore, the scalar 

principal symbol of C/g,jv is given by ■m{g) = {A4{g)eA, ca)""^ = 2~2 (a ^^-^i^j) — c)) 2 [Z, §5]. This factor is 
cancelled in the trace operation, leaving 

TrU N = - gj7r(mn+a)((m,n),(l-g)(m,Ti))_ 

The traces Trllgj^ can be analysed in the same way because in this example Ug^j^ = Ugt^pj. Indeed, 
g Ug^N is the metaplectic representation of SL(2, 2ZjNZ) realized on spaces of theta-functions [Z]. This 
connection makes it possible to get exact results on the level spacings of Ug^jss even though it is a quantum 
chaotic map. On the other hand, the results are quite different from the GOE behaviour conjectured for 
generic quantum maps (cf. [Kea]). 

Remarks on the trace expansion We do not actually use the trace expansion in this paper. This is for 
two reasons. The main one is that due to the eigenangle rescaling the significant powers U^^^ of the quantum 
map in the formula for the PCF are those for which t is on the order of IS! . But the trace expansion as it 
stands is only an asymptotic expansion as ^ 00 for fixed I. 

The'standard wisdom' regarding exponential sums (see [B.I] [II.1][G.K]) is that the van der Corput method 
produces a dual exponential sum which is only simpler than the original (in general) when the dual sum has 
fewer terms. In the case of the trace expansion, this gain in simplicity only occurs if the number of fixed 
points of is slowly growing in I. Only if the topological entropy of 4> equals zero can this be the case. The 
most favorable case is probably that of completely integrable systems, for which the £th term has roughly 
I critical points. Then the periodic orbit sum contains only I terms. But since the significant terms occur 
when I ^ N this is no simplification. 

The trace expansion is not necessary in the integrable case because the eigenangles can be written down 
explicitly by the WKB method. This is the essential use of complete integrability in this paper and also in 
[Sa.2]. The chaotic case is much more difficult because no explicit formula for the eigenangles is known and 
because the trace formula leads to a dual sum with an exponentially growing number of terms. 

3 PCF for Hamiltonians: Proof of Theorem A 

As mentioned in the introduction, the pair correlation problem for quantized Hamiltonians is similar to that 
for ZoU surfaces as presented in [U.Z]. We can therefore follow the exposition in [U.Z, §3] to the extent 
possible and omit details which are essentially similar to the Zoll case. 

There is no difference in this problem between the case of M = CP^ and the general case of symplectic 
surfaces. So in this section M can be any closed symplectic surface. However, we will make some generic 
simplifying assumptions on the Hamiltonian. The first one is 

Assumption MORSE: H : M is a Morse function. 

Let £ denote the set of values of H and let ci < C2 < . . . < cm+i denote its set of critical values. Then the 
inverse image H~^{c^, c^+i) consists of a finite number N{i') of connected components Xj each diffeomorphic 
to (cj,, Ci,+i) X S^. Hence for EG (c^, Ci^+i), H^^{E) fi consists of a periodic orbit 'IjiE) of expfEn- Its 
minimal positive period will be denoted by T^{E). 

For the sake of simplicity we will make a second assumption: 

Assumption Q: Tj' and Tj^ are independent over Q if j 7^ k. 

Under this (generic) condition, the manifold V of periodic points which arises in the calculation of p2 has 
the minimal number of components. In the general case, the formula for p2 is given by a sum over connected 
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components with dependent period functions [U.Z, Theorem 3.3]. In the simpler case (which exhibits all of 
the ideas) we will prove: 



Theorem 3.0.2 With assumptions MORSE and Q on H, the limit pair correlation function is given by 

M N(u) 



dE 



[ f{x)dp2{x) = Vm + E E E / f{kT^{E))T^{Ef 
for any f such that f e C~(]R). 

Proof: To determine the asymptotics of the sequence p"^\f) we form the generating function 

oo 
JV=1 

We wish to show that T/ is a classical Hardy-Lagrangcan distribution on S^. The asympotics of p2^\f) 
can then be determined from the singularity data of Tj. 

To show that ^ f{9) is a Lagrangean distribution, we write it as the trace of a Toeplitz type Fourier 
Integral operator, defined as follows: We form the product manifold X x X and consider the product Szego 
projector 

n O n : L^{X xX)^ ij2(X) O H'^{X) 

and the diagonal projector 

oo oo 

^diaa := Hjv n^v : L^iX x X) ^ Ojv Qn- 

JV=1 N=l 

We next observe that 

JTR 

Noting that N is the eigenvalue of the number operator M we can rewrite this in the form 

Trli-N^^N I f{t)e^*'^^"''^'-'^^"''Ut. 
Jm. 

The generating function is then given by 

jv=i 

= TrUMage''^^^'^ I f{t)e''^^"^'-'^^"Ut. 
Jm, 

We recall here that H = ITiJII where H is the pull back to X of the function so denoted on M. We now 
have to analyse each operator which occurs under the trace sign. 

(a) Hdiag- This operator is the composition of the product Szego projector 11 (g) 11 with the full diagonal 
weight projection 

Pdiag -.L^X^X)^ L%{X) L%{X) 
where L%{X) is the eigenspace of of eigenvalue N. 



13 



From [G.S.2] [U.Z] it follows that Pdiag is a Fourier Integral operator in the class I^{X^'^^ x X^^^ , f) where 
X(^) = X X X and where F is the flow-out of the coisotropic cone 

e = {(Ci,C2) e T*(X(2)) : //(Ci) - H{C2) = 0}. 

That is, let = exptEH x exp-tEn denote the Hamilton flow generated on T*(X(^)) by -ff(Ci) - -f^(C2)- 
Then in a well-known way, the map 

ie :IRx e^r*(X(2)) X T*(X(2)), (f,Ci,C2) - ((Ci,C2),$*(Ci,C2)) 

defines a Lagrange immersion with image equal (by definition) to the fiow out F. 

On the other hand H 11 is the exterior tensor product of two Toeplitz (hence Hermite type Fourier 
Integral) operators. According to [B.G, Theorem 9.3], we therefore have 

Ili»U = a + (3 

with 

Q;g/o(x(2) xX(2),SxS) 

and with WF{P) contained in a small conic neighborhood CofSxOUOxS. Moreover, the symbol of 11 (g) Ft 
is given by 

cr(n® n) = cr(n) ® cr(n) 

on S x S — C Hence H (g) H is essentially a Toeplitz structure on the symplectic cone S x S C T*{X^'^^ ). 
The complication due to C will ultimately prove to be irrelevant in the analysis of the trace since it will not 
contribute to the singularities along the diagonal. Hence we can (and will) pretend that this component of 
WF{Il (g) n) does not occur. 

By the composition theorem for Hermite and ordinary Fourier Integral operators [B.G, Theorem 7.5] it 
follows that (modulo the term Pdiag ° /?) 

Hd.ag er(X(2) xX(2),F) 

where F is the fiowout Lagrangean in S x S for the co-isotropic subcone := fl S C S. That is, the map 

ie ■■= «ekxe : IR x 6 ^ S x S 

is a Lagrange immersion with respect to the symplectic cone E x E and F is its image; of course it is only 
an isotropic immersion with respect to T*{X^'^^ x X'^^). 

(b) ei«[-^®^l: This operator does not require a fancy analysis since Af is simply the differentiation operator 
by the generator ^ of the contact fiow. Hence e*^!-^®^' is the translation operator F{x, y) F{(1>^ {x),y) by 
(j)^ xid onX(2). 

(c) H (g) H J^ f{t)e'^'-^"'^'-''^^"^dt : Here we have inserted the factor H O H, as we may, to simplify the 
discussion. 

Since [Af, H] = = [Af, H], we may remove the projection H from the exponent. Then the unitary under 
the integral is given by 

Each factor is the exponential of a pseudo differential operator of real principal type and is therefore Fourier 
Integral. It follows again from the composition theorem [B.G, Theorem 7.5] that 

jjgitffAT g /-i (X X X, (7 n H X S X S), C :={{t,T,x,^,y,v):T + <jm{x, OH{x) = 0, ^\x, = {y, v)} 
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where tp* is the Hamihon flow on T*X generated by a_s/{x, S,)H. Note that aj^^^^^^ = -j^), which generates 
the hft of the central circle action on X to T*X . Also, the Hamilton flow of H on T*X is the two- fold lift 
of the Hamilton flow of H on M: first from M to X and then from X to T*X. Since the Hamilton vector 
field of aj^(^x,i)H on T*X is given by 

and the Lie bracket of the two terms is zero, we have 

V-* = expfiJS<^j^(^ ^) oexpt(Tjv-(a;,OSjf 
The isotropic cone C fl S x S can be parametrized by 

ic^ : M X E ^ IR X E X S, (t, C) ^ {t, a^HiQX, V'*(C))- 
The symbol of H O H /^^ f(^^^git(AfH0i-i0AfH)fj^^ niay then be identified with the spinor f{t)\dt\ ^ (gian- 

Putting the above together we see that up to the factor of Pdiag the operator under the trace is a Hermite 
Fourier Integral operator associated to the graph of (j)^ x ido x 'tp~* on S x S. The effect of the Pdiag 
factor is to reduce this torus action to the quotient of S x E by the diagonal contact flow. 

The remainder of the calculation is similar to the ZoU case, so we will just summarize the key ideas and 
refer to [U.Z] for details. First, the singularities of T are caused by the periodic orbits of the reduced flow 
exptEn X exp—tEn on the characteristic variety V := {H{zi) — H{z2) = 0} C M x M. The singularities 
occur at the holonomy angles of the lifts of these periodic orbits to the prequantum bundle of M x M. 
Since the periodic orbits of the product (or difference) flow are products of periodic orbits of the factors, 
the holonomy of the periodic orbits are ratios of holonomies of the factors. In particular, if for each period 
there is just one periodic orbit, then the ratio is always one and there is only a single singularity at ^ = 0. 
The general case is discussed in detail in the analogous situation of [U.Z]. 

The set of relevant periodic points and their periods is given by 

V = {{zi,Z2, T) e V X suppf : expTSH(zj) = zj}. 

Observe that each component of V is of dimension 3: the T = component equals V, and the others consist 
of continuous unions of products a x (3 oi periodic orbits a, (3 with the same (or rationally related) periods 
(see [U.Z, Lemma 3.7]). Under assumption Q, only 'squares' a x a of periodic orbits arise. It follows that 
each component contributes with equal strength to the pair correlation function. 

To complete the calculation we need to determine the canonical density on V whose integral gives the 
principal symbol of T j at = (or at other possible singular angles) . On the T = component it is 
just the Liouville density of V and hence the contribution of this component is the volume of V. For the 
T > components we use the Morse theory of H to break up M into local action-angle charts XJ. We 
then wish to determine the surface measures on the surfaces V fl {Ty{z\) = Ty{z2)} C x X^ where as 
above T^{z) is the minimal period of the periodic orbit thru z. We will express the result in terms of the 
local action angle coordinates {Ii, (pi, I2, (P2) on Xj x Xj. The coordinates {Ii, (pi, (^2) are independent on 
V n {Tj{zi) = Tj{z2)} and (dropping the sub and superscripts) the Liouville measure in these coordinates 
equals T{dli A d(t>i A d(t)2) [U.Z, Lemma 3.11]. On the level sets T{zi) = T{z2) the sm-face measure is then 
given by T"^ d(pid(p2) ■ The measure of the T > component of V is then the integral over the set 
of periods T of this form times dT. Changing variables from T{E) to E, summing over the components X^ 
(and reinstating the sub and superscripts) gives the stated formula. ■ 
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4 PCF for quantized perfect Hamiltonian flows on CP^: Proof of 
Theorem B 



In this section we consider general Hamiltonians H on CP^ which are perfect Morse functions. Our purpose 
is to show that the pair correlation functions of their quantized Hamiltonian flows are Poisson on average. 

4.1 Classical Hamiltonian 5*^ actions and perfect Morse Hamiltonians on (DP^ 

Up to symplectic equivalence, CP^ carries a unique Hamiltonian action. The model example is that of 
rotations of 5^ C around the cca-axis, whose Hamiltonian is the coordinate. 

In general, a function I : <CP^ — > IR is called an 'action' variable if its Hamilton flow exptS/ is 2it- 
periodic. Any such global action variable on (EP^ must have precisely two critical points which arc fixed 
points of its flow. This flow has a global transversal connecting the flxed points. The travel time from this 
transversal defines an angle variable 9 symplcctically dual to I and the transformation x(^, 0) = {xz,0) is a 
global symplectic transformation to the model case. 

Now suppose that H : CP^ ^ IR is a perfect Morse function, with a non-degenerate minimum value 
equal to —1, and a non-degenerate maximum value equal to 1. Then let denote the distribution function of 
H : CP^ — > [—1, 1], i.e. ijl{E) = \{H < E}\ where | • | denotes the symplectic area. Under our assumptions it 
is a strictly increasing smooth fimction on (—1, 1). Let / : (CP^ [0, Att] be defined by I{z) ~ ii{H{z)). Then 
/ is an action variable for H: as functions on the symplectic CP^, / and H Poisson commute, {I,H} = 0, 
and the Hamilton flow of I is 27r- periodic. It is obvious that / generates the algebra a = {/ G C°°(CP^) : 
{/, H} = 0} and hence we may write H = where <f) is the inverse function to /i. Non-degeneracy of the 
critical points forces d(j) ^ since 

d^H = (j)'{l)(fl + cj)"{l)dl(^dl (13) 
so that d^H = (j)"{l)dl ® dl at critical points. 

4.2 Quantum actions over (DF^ 

It is a general fact that compact Hamilton group actions can be quantized [B.G]. In the model case, the 
quantization is given by the generator of a linear action, such as rotations around the .T3-axis (if we think 
of CP^ as the embedded S"^) . In general (without conjugation to the model case), it can be constructed as 
follows: The Toeplitz quantization H/H of 7 is a positive operator satisfying e^'^'-'^P^^l = I + K where K is 
a Toeplitz operator of order —1. In a well-known way [CV], we may add lower order terms to H/H to arrive 
at a positive operator / satisfying [/, A/] =0 and e^'^^^^ = I. Therefore Nl has only integral eigenvalues. 
Since AA = iV on H'^{N), it follows that 

Sp{i\H^iN)) = {j^:j = -N,..., N}. (14) 

4.3 WKB for quantum completely integral Hamiltonians in genus zero 

We now come to the deflnition of the quantum systems which concern us. In the homogeneous setting we 

say: 

Deflnition 4.3.1 A quantum completely integrable Hamiltonian is a Toeplitz Hamiltonian H on H^{X) 
which is given by a poly homogeneous function 

H = $(/) ~ (/.(/) -(- </._i(/)AA-i ^ (j>_2{i)N-^ + . . . 
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of a global action operator i , i.e. a generator of a quantum Hamiltonian action. Automatically, [H,J\f] =0 
and we write H^^^ for UnHUn- 

The principal symbol defines a perfect Morse function on <CP^ if (j)' ^ 0. Conversely, suppose that H 
is a perfect Morse function on M, so that there exists a global action I and monotone (j) for which H = 4>{I). 
The Toeplitz quantization litqHUjq of H does not necessarily commute with the quantization I of /, so it 
is not quantum completely integrable in the strong sense of the definition above. To make it commute it 
is only necessary to average Ylj^HYli^ with respect to the quantum S'-'^-action expitl. Since this docs not 
change the principal symbol it may be regarded as a quantization of H and we will reserve the notation H 
for it. Since H commutes with I and since I has a simple spectrum, H = $(7) for some function It 
must be a polyhomogenous function since if is a Toeplitz operator. We therefore come back to the same 
definition as above. 

The following proposition is an immediate consequence of the polyhomogcneity: 
Proposition 4.3.2 Denote the eigenvalues of H'^'^^ by Sp{H'^'^^) = {Xnj}- Then: 

4.4 The quantizations U-^t,N and Ut,N 

There arc two approaches to the quantization of a Hamiltonian flow xi = cxp t'E.H '■ (i) by first exponentiating 
and then quantizing or (ii) by first quantizing and then exponentiating. As above we would like to produce 
a quantum completely integrable system, so that U^^ should commute with I. This could be done by 
quantizing H ^ H as above and then exponentiating to get C/t,jv = Hn^^^'^^^n- This automatically 
produces a quantum completely integrable unitary group. 

Alternatively, we could exponentiate H to get the flow exp tSjj and then quantize the flow by the Toeplitz 
method. By Proposition 2.3, exptSn lifts to a contact transformation xt on X so by [Z, §3] there exists a 
canonical symbol at S S'^ {T* (X)) such that the L^^t.^ ^N<^tXt^N is unitary. There is no guarantee that 
the resulting quantization commutes with the quantum action generated by / nor that it is a unitary 
group. The first defect can be overcome, as above, by first averaging against the action and then applying 
the functional calculus as in [Z] to make the quantization unitary. We henceforth use the notation C/^t.JV to 
refer to the result. 

To compare the methods wc prove: 

Proposition 4.4.1 Ut,N — Uxt,N modulo operators of order —1. 

Proof: To be cautious, we first convince ourselves that IItv and e''*-'^^^^ commute. Indeed, IIjv = 11 o 
Jgi Ad{e^^^)e~^'^^ dO and the averaging operator commutes with both 11 and multiplication by H. It follows 
that 

TT* TT — TT „-itJ\fnHnTT AtMnHnTr _ TT „-iWnHnAtj\fnHnTT _ tt 

Thus, Ut,N is unitary on H'^{N). 

Wc then observe that Ut = He^^-^^^^Tl = Ile^'-^^II is a unitary group of Fourier- Toeplitz integral 
operators whose underlying canonical relation is the lift of the Hamilton flow of H to the symplectic cone 
S C T*{X) — {(x,^,ra) : r G M^}. This follows very much as in the proof that exponentials of first order 
pseudodifferential operators are groups of Fourier integral operators. Indeed, IU\fHIl is a first order Toeplitz 
(pseudodifferential) operator with principal symbol rH. Hence its exponential is a Toeplitz Fourier integral 
operator with bicharacteristic fiow equal to the Hamilton flow of rH on E. By Proposition 2.3, this Hamilton 
flow is the lift of exptSij to S. 
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On the semi-classical (non-homogeneous) level, we note that Ut^N = Uns^^'^^^^^n forms a non- 
homogeneous Hermite- Fourier integral distribution as N varies. By non-homogeneous is meant an oscillatory 
integral (with complex phase) associated to a non-homogeneous Lagrangean; here it equals the graph of Xt 
on X. Both Ut,N and Uxt,N are Hermite-Fourier integral distributions associated to the same Lagrangean, 
the graph of Xt^ and both have principal symbols equal to the graph 1/2-density. It follows that they can 
only differ by operators in the same class and of order -1. In fact both are functions of / of the form 
exp(i$t(j). The only possible difference is in the subprincipal terms in the polyhomogeneous expansions of 
the $ functions. ■ 

It seems most natural to concentrate on the unitary groups Ut,N so we will always assume our quantum 
maps have the form expjt<E>(/) for some polyhomogeneous It is straightforward to extend them to the 
more general quantum maps U^^^n- 

4.5 Final form of TrC//jv 

It follows from the above that the pair correlation function for the completely integrable quantum map 
exp{itJ\f (f)(1)) is given by 



pifif) = m + ^T.h^^)\stiN,i)\' 



where 7^ and where 



N 

5t(A^,£) =^e*"^*(*). (15) 



We will refer to this as the 'homogeneous' case since the Hamiltonian has the form .H'(^) = 7V0(/(^)). 
In the general polyhomogeneous case, the exponential sum has the form 

N 

St(7V,£) = ^e**^^*(*'^) (16) 
i=i 

with <I>(-i-,A^) = (/>(-^) + + .... For the sake of simplicity wc often assume the Hamiltonian is 

homogeneous, but we would like to note that they also hold in the general case since the lower order terms 
do not affect the lattice point counting arguments that are the heart of the matter. 

As regards the individual sums we also note that the terms of order or lower in $(j, N) cannot affect 
the pair correlation function. That is, let us put: 

Z,(iV,^) =f;e^*^^*=(^'^) (17) 

with $2(i,A^) := [0(^) + ^^77^ + ^^^]- The following is obvious: 
Proposition 4.5.1 For f with f G C~(IR), 
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4.6 Proof of Theorem B: pair correlation on average 

We now come to the main results. In the first we allow for general polyhomogeneous exponents. 

Theorem 4.6.1 Let Ut,N = exp{itAf^{I)) as above, with \4>'{x)\ > Ci. Then the limit pair correlation 
function p2,t for Ut,N is Poisson on average.- 



^£p2,t(/)di = /(0) + /(0) 



for any interval [a, 6] C IR and any f G . 
Proof: By the above, it suffices to show that 



N „ ,„ 

To prove this, we use the Hilbert inequality (cf. [Mo, §7.6 (28)]) 

.b N N 



^ .b N N 

±_ / |^e2"''^*|2dt = iV + 0(^^) 

' ^■^'^ .7=1 .7=1 '^^^ 



with 



5j = min - 



The two terms correspond respectively to the diagonal and to the off-diagonal in the square, and the 0-symbol 
is an abi 
so that 



is an absolute constant (which can be taken to be 3/2). In the case at hand, = Ni[(j){j^)-\- ''' '^""^ -t- *^ 



5N,j > e [min^e[a,b0i^)\ + 0{1/N)]. 
It follows that i{\(f)'{x)\>C>0 then 



1 N 

\St{N,e)\^dt = N + 0{-). 



b — a J„ £ 

Therefore the limit equals 



lim — > /(——)+ 0(limsup — > -rfi^rr))- 
The first term tends to /j^ f{x)dx = /(O). When supp / C [— C, C] the second is 



« ^[ E i\ « 



logN 



e<CN 
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4.7 Proof of Theorem B (a): Number Vciriance on average 

Next, we prove that the number variance is Poisson on average: 
Theorem 4.7.1 With the same hypotheses as above, we have 

rb 



hm / Y.i^,ht = L. 

N^oo b-aj^ 2,t 



Proof 

We have 



i=l e=l j=Ak,i,k=l 



The first term is a Riemann sum and we have 



. ^,Tr£L, f .smLx.n , 

( ) ax = L. 

X 



The second can be expHcitly evaluated as above and is bounded by 



P N ^ ^ N 



5 Mean square Poisson statistics for quantum spin evolutions: 
Proof of Theorem B (b) 

We have just seen that averages of the pair correlation function pj^-* of quantized Hamilton flows Ut^N 
converge to pPOissoN ^ loixg as -ff is a perfect Morse function. In particular, the result is true for linear 
functions H = al in the action. Since exponential sums with linear phases are far from random (see §6), it 
is evident that the averaging is the agent producing the random number behaviour. 
A much stronger test of the Poisson behaviour of Ut^N is whether the variance 



^l\pZ\f)-P2'''''''''{mt 



tends to zero as N ^ oo. It is easy to see that quantum Hamiltonian flows with linear Hamiltonians do not 
have this property (§4.1). We therefore turn to quadratic Hamiltonians H = + f3I. Our next result shows 
that the pair correlation functions of their quantized Hamilton flows converge in mean square to Poisson. 
The proof is based on techniques from mean value estimates on exponential sums; see [B.I] [G.K] [H.l] for 
background. 

Theorem 5.0.2 Let Ha,i3 = oil^ + 1^1 and let p'^f^^ ^ be the pair correlation measure for the quantum map 
U{t,a,i3),N = ^xp{itJ\fH(^a,i3;N))- Then for any t ^0, any T > and any f G <S(IR) with f G C^(IR) we have 



)• 
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Proof: Removing the £ = and diagonal terms as above, it suffices to show that 

llfjw^fiw^l t ««(«^^ + W-*)))P"^-0(<!5I^)^ (20) 

To prove this, we use the BeurUng - Selberg function S^jji^. It has the following properties (see [G.K]): 

• BT,\t\5 > X[-T,T]'i 

• Sn^vBT,\t\sC{-\t\5,\t\5). 

Here, X[-t,t] is the characteristic function of [— T, T]. Then the integral on the right side above is bounded 
above by 

(■ r ^ ^ f JL •7-2 _ p 

/ / BrMs{a)Br,it\sm-^J2f(N^^ E e{t{a^——+P{j-k)))fdadp. (21) 
Squaring and evaluating the Fourier transforms gives 

1 0—1) ^ ^ '1 1 1 -I 1 1 

]^ E E /(|)/(|)[ E E S^,W.(i(^l^^-^2^^))BT,|*|.(i(4(il-fcl)-^2(i2-fc2))). 

^1#0£25^0 jlT^fcl Jl,fel = l j25^fe2 J2,fe2 = l 

(22) 

By the support properties of -B^ the latter expression is bounded above by 

]^EEi/(|)ii/(|)i^w^i'^^) (23) 

with 

•2 iJl '2 7 2 

I{N,h,i2) = #{{ji,ki,j2,k2) G [l,N]^:ji^ki : \{£,ilj^-i2^-2-Z-^)\ < S,\{£,{j,-k,)-£2{j2-k2m < S}. 

(24) 

Introduce new variables hi = ji — ki,mi = ji + ki so that the conditions read 

\£i^-£2^)\<5 

\£ihi-£2h2\ <5 

The change of variables is invertible so /(iV, £1,^2) is the number of integer solutions (/ii, /12, Wi, TO2) with 
\mi\ <2N,\hi\ <iV- |mi|. 

Since \£ihi — ^2/12] € N it can only be < 5 if it vanishes. Therefore the second condition is equivalent to 

{£,h) = ^h2 = ^^hi. 

t2 

Here wc assume \hi\ > |/i2| so that = max{\hi\, |/i2|} ^ \h\ and we abbreviate £ = (^1,^2) etc. Substi- 
tuting in the first condition we get 

£i/ii(mi -TO2) =0{N6). 

Now let us count solutions. We split them up into two classes: (a) homogeneous solutions with mi — 
7712 = and (b) inhomogeneous solutions with mi — m2 7^ 0. The homogeneous solutions are lattice points 
{hi, mi,£i;h2, m2, £2) which solve the system 

£ihi=£2h2 (/ii^0,/i2^0); 
mi = m2. 
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Clearly there are 2N solutions of mi = ni2- For each integer s in [—N"^, N"^] there are < d{s) ways of writing s 
as a product hi£i with hi,£i G [—N, N]. Here, d{s) is the divisor function (the number of non-trivial divisors 
of s). Hence the number of solutions of the homogeneous system is 0{N J2^=i ^(■s)^) = 0{N^{logN)^). 

Now let us count inhomogeneous solutions. We write Iih,{N,ii,£2) for the number of inhomogeneous 
solutions with fixed (^1,^2)- Since /12 is determined from (ft,i,mi,m2) it suffices to count these triples. First, 
there are 0{N) choices of mi. Then put m2=mi+M with M > 1 so that hiM = 0{^S). From M <2N 
the number of pairs {hi , M) is bounded above by 

Hence/,,, (TV, £1,4) « 1^1^. 
It follows that for / e C~(IR), 

^Y.e.^,Y.e.^,\K%)\mmNJiM<<0{^^) + '^^^^^ 

(25) 

« 0( (1^^) + o{^^) = o(M!)). 



5.1 Proof of Theorem B(b) for general non-degenerate phases 

We now consider a much more general class of Hamiltonians for which a similar result holds: 

Theorem 5.1.1 Let 13 = oi<j){i) + (31 where > C > on [—1,1] and let P^^(^^ ^ 

correlation measure for the quantum map U(t,a,p),N = ^xp{itMH{^a,i3;N))- Then for any t ^ 0, any T > 

and any f e <S(]R) with f e C^(]R) we have 



Proof The previous argument now leads to the lattice point problem: 

ilhi = i2h2 

him) - <t>i^)) - ^-^m) - H^)) = o{i/N) 

Here as above hi = ji — ki. By the mean value theorem there exist ^j^kt € [ji, fci] such that 4'{^) — 4'ij^) = 
~ ^i)- above we then get the system of constraints: 

hhi = £2h2 

hi£i{cl>'i^j,kJN) - cj^'i^i^kJN)) = 0(1) 
Then writing (/>'(^iifci/^) - ^'(Osfcs/^) = 0"(^iifcii2fe2/^)(^jifci - ihk^)/^ we get the system: 

iihi = £2h2 

hih<i>"{^hk.nkJN){^j,k, - ^j.k.)) = 0{N) 
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By the assumption |0"| > c > this gives 

iihi = e2h2 (ej^ojij^o) 

(26) 

Mi(e,ifci-e,.fej = o(iv) 

Let us change to the variables {hi, rrii) as in the quadratic case and write ^j.fe. — ^{hi, rrii). We wish to count 
the number of lattice points hi,mi, £2, ^2,1112) satisfying the system 

iihi^i2h2 iij^O,hj^O) 

(27) 

\ahi,rm)^ah2,m2)\ < 

As in the quadratic case, we regard (/ii,^i,mi) as independent variables, so that the first equation is a 
constraint on (/i2,^2)- We also regard the second constraint |^(/ii,toi) — £,{h2,m2)\ < as a constraint 
on \m2\. To put it in a more convenient form we consider ^(/i2,'Ti2) as a function .^^.2(7712) and invert the 
function ■ An easy calculation gives 

^,(m) = Ni^TH f + (2^ - 1)^)^^) (28) 



hence 



A^^(^) ^ C <t>'{— + (2s - l) — )ds)}\ r (/."(— + (2s - l) — )ds]. 

dni ' ' ^ ^Vo 27V ^ '2N' ^^7o 27V ^ '2N' ^ 



(29) 



Since C < \4>"{x)\ < C" for x e [—1, 1] and certain positive constants C, C", It follows that that |^^/i(™)| > 
(5 > (Vm G [0, 2A^]). Therefore a smooth inverse function ^^^^ exists on the range of ^. In particular, 
is Lipshitz, so (|27|) is equivalent to the system 



£1/11 =£2/12 (^, ^0,/i, ^0) 
|m2-C/(e(/ii,mi))|<C"^. 



(30) 



The situation is now very close to the quadratic case: There are d(hi£i) solutions (/12, ^2) of the first equation 
and N values of mi. For given (ft-i, ^i, /12, ^2, wi) there are < 1 + 0{-fp^) solutions m2 of the second 
equation. Summing 1 over the relevant {hi,£i, /i2,^2,™i) gives 0(A^^(log A^)^ as in the homogeneous part of 
the quadratic case. Summing 0{j^f^ over these lattice points gives 0(A^^(log7V)^) as in the inhomogeneous 
part of the quadratic case. ■ 

Corollary 5.1.2 Let Nm = [m(log m)**] ([■] = integer part). Then for almost all {a, (3) with respect to 
Lebseque measure and all t ^ we have 

Proof By the above, 

^T^vl I ipSi.)(/)-pr^^"""(/)p^«^/3<oo. 

m=l ^ ' J-T J-T 

Since the terms are positive it follows that for almost all (a, /3), 

E IP^S.,)(/) - P2^°^""°^(/)Pd«rf/3 < 00 

m—1 
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and for these {a, (3) the mth term tends to zero. The set of such {a, (3) apriori depends on the smooth 
function /. However by a standard diagonal argument one can find a set of full measure that works for every 
continuous / (see [Sa.2][R.S] for further details). 

Remark In this corollary we have adapted an argument from [Sa.2][R.S], where the pair correlation problem 
is studied for flat tori and for some homogeneous integrable systems. Their main result was that the relevant 

pair correlation functions are almost everywhere Poisson. After proving the almost everywhere convergence to 
Poisson along a slightly sparse subsequence (as in the above Corollary), they show that for iV^ < M < Nm+i, 



as m — > oo. This last step seems to be much more difficult in our problem. The difference is that the 
spectra in [Sa.2][R.S] increase with increasing A'' and the common terms cancel in the difference above. On 
the other hand, our spectra change rapidly with N and there are no (obvious) common terms to cancel. In 
[Z. Addendum] we will show that an a.e. result comparable to that of [R.S] holds for the average in N of 



6 Appendix: Linear and quadratic cases 

In the case of linear and pure quadratic Hamiltonians, the exponential sums discussed above are classical and 
there are many prior results in the literature. We briefly discuss what is known and add a few observations 
of our own. 

First, the pair correlation problem for linear Hamiltonians H = al has been studied since the fifties. Sec; 
[B1.2][R.S][Ca.Gu.Iz] for discussion and references to the literature. The main result is that only three level 
spacings can occur for a given a and the pair correlation function is not even mean square Poisson. 

In the case of quadratic Hamiltonians, we get the incomplete Gauss sums: 




(AT) 
^2;(t,a,/3) 



if the phase is a polynomial. 



JV 



St{N;£) = ^e'^^itNe[(jrf+ajr]_ 



i=i 



In the special case a = and t = 1 they are classical complete Gauss sums 



N 




If {£, N) = l then 



\G{i,0,N)\ = 




l{mod2) 
0(rnorf4) 
2(mod4) 



In general 



G{i,0,N) = {i,N)G{ 



)• 



Hence the values of 




depend on the residue class of N modulo 4. If A' = 2 (mod 4) then 1^ = 0. If A' is odd, then 





ejta 
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When N = p, a, prime number. Then {i,p) = 1 except for multiples kp with k G supp{f). They make a 
neglible contribution, so 

P P Jm 



The sums 



Thus the prime sequence is Poisson. 

At the opposite extreme, suppose A'' = p'^ for some prime p. Then 

^ k=0 q=l,(q,p) = l ^ 

1 ^""^ ^ g 

pM-k E f^pM-k ) 

9=1,('2,P) = 1 

have the form of Riemann sums for /(O) as M — fc ^ oo except that the terms with p\q are omitted. These 
terms also resemble Riemann sums for /(O), multiplied by 1/p. Hence each term is roughly 1 — 1/p times 
/(O). Since there are M such terms, the coefficient of /(O) tends to infinity and the pair correlation function 
cannot be Poisson. 

6.1 Poisson on average in t 

If we allow t to vary then we do have an average Poisson behaviour: 

Proposition 6.1.1 For any interval [— T, T], the average PCF of IlexpitJ\f P is Poisson, i.e. 

7V2 



Proof: The integral equals 



N 



N ^ e-'^N' ^ f-k 



2 



Z^J^n'^Z^ Nhmi 

ijiO TO=1 0<|/t|<iV-|m| 

where as above h = j — k,m = j + k. Using just that sin a; << 1 this is 

]vE-^'(]v^ E E 

£^0 m=10<\h{<N-\ml 



We have not determined whether the variance tends to zero in this case. 
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